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Abstract 



Semiclassical Einstein-Langevin equations for arbitrary small metric perturba- 
tions conformally coupled to a massless quantum scalar field in a spatially flat 
cosmological background are derived. Use is made of the fact that for this prob- 
lem the in-in or closed time path effective action is simply related to the Feynman 
and Vernon influence functional which describes the effect of the "environment" , 
the quantum field which is coarse grained here, on the "system", the gravita- 
tional field which is the field of interest. This leads to identify the dissipation and 
noise kernels in the in-in effective action, and to derive a fluctuation-dissipation 
relation. A tensorial Gaussian stochastic source which couples to the Weyl tensor 
of the spacetime metric is seen to modify the usual semiclassical equations which 
can be viewed now as mean field equations. As a simple application we derive 
the correlation functions of the stochastic metric fluctuations produced in a flat 
spacetime with small metric perturbations due to the quantum fluctuations of 
the matter field coupled to these perturbations. 
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I. INTRODUCTION 



The semiclassical equations of gravity are the generahzation of Einstein equations when 
the matter source is described quantically. The source of these equations is the vacuum 
expectation value of the stress tensor of the matter fields. This theory which assumes that 
the gravitational field is classical should work well outside the Planck scales and when the 
quantum fluctuations are small because the matter source is an average term. But as it has 
been emphasized by several authors one should expect that, if the fluctuations above 
the average are large, and this depends on the quantum state of the fields, the semiclas- 
sical equations should not give a correct description. In fact, one expects that a better 
approximation should describe the gravitational field in a probabilistic way. In other words, 
the semiclassical equations should be substituted by some Langevin type equation with a 
stochastic source which describes classically the quantum fluctuations. A significant step 
in this direction was made by Hu 0] who proposed to view the semiclassical back-reaction 
problem in the framework of quantum open systems, where the quantum fields are seen 
as the "environment" and the gravitational field is seen as the "system". Following this 
proposal a systematic study of the connection between quantum field theory and statisti- 
cal mechanics has resulted in the development of a consistent mathematical framework in 
which the relevant semiclassical Einstein-Langevin equations can be derived @] The 
key to these results is the influence functional method used in nonequilibrium statistical 
mechanics to describe the environment-system interaction when only the state of the system 
is of interest. Here, the environment are the quantum fields whose degrees of freedom are 
not of interest in the semiclassical back-reaction problem and are traced out, i.e. coarse 
grained, and the system is the gravitational field which is the field of interest. The infiuence 
functional provides information about the dissipation and the noise suffered by the system, 
and the corresponding fluctuation-dissipation relation. In this framework it is seen that 
external stochastic sources with a given probability distribution appear in the semiclassical 
equation for the gravitational fleld. The origin of the noise acting on the system are the 
coarse grained quantum flelds whose degrees of freedom have been traced out. 

In this paper we derive the semiclassical Einstein-Langevin equations for arbitrary linear 
metric perturbations conformally coupled to a massless quantum scalar fleld in a spatially flat 
cosmological model. Semiclassical Einstein-Langevin equations have been recently derived 
by Hu and Sinha @] for small anisotropics conformally coupled to massless flelds in a spatially 
homogeneous background working in the framework of quantum cosmology, and by Hu and 
Matactz |]5| who derived the semiclassical equations for the scale factor in a spatially flat 
universe due to the coupling of different quantum scalar flelds. 

We work in the "in-in" or CTP (closed time path) effective action framework. The CTP 
effective action technique which was flrst proposed by Schwinger is an effective action 
method adapted to compute expectation values of quantum operators rather than matrix 
elements as in the usual "in-out" effective action method. It has been adapted by Jordan 
and Calzetta and Hu to a curved background and from it real evolution equations which 
admit an initial value formulation can be derived. It has been used to derive the semiclassical 
equations for small anisotropics in an homogeneous background due to quantum matter flelds 
0] [1^ and it was also used in a previous paper |Tl| , referred from now on as Paper I, to derive 
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the semiclassical equations for arbitrary linear perturbations in homogeneous backgrounds 
due to massless conformally coupled matter fields. In this last case the CTP effective action 
method leads in a very direct way to the evaluation of the vacuum expectation value of 
the field stress tensor which had been obtained previously by other methods The 
CTP effective action method is also currently used to study dissipative effects produced by 
the inflaton oscillations around a true vacuum []T3[. These effects are of great interest to 
understand the reheating mechanism in the inflationary cosmological scenario. 

The connection between the CTP effective action and the influence functional was es- 
tablished by Calzetta and Hu |^ in the semiclassical context. It turns out that in this case 
there is a very simple linear relation between the so called influence action and the CTP 
effective action. That relation is not so simple if one quantizes the gravitational fleld or one 
wants the fleld equations for the mean quantum fleld on a given background |14] since then 
one has to advocate decoherence to go from a quantum to a classical description. This, of 
course, is of great interest for the study of structure formation in the early universe. In our 
semiclassical case, however, the CTP effective action already partially derived in Paper I 
gives all the necessary information not only on the dissipation effects but also on the noise 
due to quantum fluctuations. All we need to do is to identify the corresponding dissipation 
and noise kernels. 

The plan of the paper is the following. In Sec. H we briefly review the CTP functional 
formalism and derive the complete CTP effective action for arbitrary linear metric pertur- 
bations in a homogeneous background. Note that in Paper I the CTP effective action was 
only partially given since the complete action was not necessary to derive the (mean fleld) 
semiclassical fleld equations. In Sec. |IT1| the influence action is derived, the dissipation 
and noise kernels are obtained and the fluctuation-dissipation relation is given. In Sec. ^ 
we derive the semiclassical Einstein-Langevin equation and identify the stochastic source. 
We show that a tensorial Gaussian stochastic source which couples to the Weyl tensor of 
the spacetime metric is the origin of noise in this case. Finally, in Sec. |V| we apply the 
Einstein-Langevin equation in a non cosmological background and derive the correlation 
functions of the stochastic metric fluctuations produced in a flat spacetime with small met- 
ric perturbations as a consequence of the coupling of the quantum matter fleld to these 
perturbations. 



II. CTP EFFECTIVE ACTION FOR WEAK INHOMOGENEITIES 



First, let us sketch the CTP functional formalism for the evaluation of the CTP effective 



action; for a more detailed exposition see Refs. p |n[] [jT5[. Let us consider the quan- 
tization of a scalar fleld 4>{x), the usual in-out effective action is based in the generating 
functional W[J] which is related to the vacuum persistence amplitude in the presence of some 
classical source J(x) 
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From this functional one generates matrix elements of the fleld (f){x) 
by functional derivation with respect to J{x). Now in order to work with expectation values 
rather than matrix elements one deflnes a new generating functional whose dynamics is de- 
termined by two different external classical sources J+ and J_ by letting the in vacuum evolve 
independently under these sources, i.e. exp{iW[J+, J-]} = X]a(0, T) j_ (a, T|m, 0) j^. 
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where we have assumed that {\a, T)} is a complete basis of eigenstates of the field operator 
4>{x) at some future time T. This generating functional admits a path integral representation 
which we may write in a compact form as 

^iW[j+,j.] ^ f j^^^^^j^^^_y{S[M'Si^-]+J+<l>+-J-'P-)^ (2.1) 



where J±4>± stands for / d"'xJ±{x)(j)±{x) and it is understood that we sum over all fields 0+ 
and (j)- with negative and positive frequency modes, respectively, in the remote past (these 
are the in boundary conditions) but which coincide at time t = T (in practice one usually 
takes T —>■ oo, the remote future). Thus this path integral can be thought of as the path 
sum of two different fields evolving in two different time branches: one going forward in time 
in the presence of J+ from the in vacuum to a time T, and the other backward in time in the 
presence of J_ from time T to the in vacuum, with the constraint that 0+(r, x) = 0-(T, x). 
For that reason this formalism is called CTP functional formalism. Note that we use an 
arbitrary number n of spacetime dimensions in order to perform dimensional regularization. 
From this generating functional expectation values can be obtained. Let us define 

^'^l^-'-J.±fclJ,.J_l. (2.2) 



6J. 



± 



and assume that these equations can be reversed, then the CTP effective action is the 
Legendre transform of the above functional 



CTP 



W[J+, J_] - J+0+ + J_0_, (2.3) 



where it is understood that the external sources are functionals of the fields 0+ and 0_ 
through the definitions ( p^.2[ ). From this by functional derivation with respect to 0± we get 
the equations for the expectation values 4>±[J+, J-]- The equation for the vacuum expecta- 
tion value of the field 0[O] = 0±[O,O] = (0, m|0(x)|m, 0) is then simply obtained imposing 
that J±. = 0: 



6T 



CTP [<P+, ' 



= 0. (2.4) 

4>±=m 

To simplify the notation it is useful at this stage to introduce a more compact notation 
<S[4>a] = S[(j)-^] — S[(l)-], (/)a{x) = (0+, (j)-), and Ja{x) = ( J+, J„), where a and b take the two 
values + and — , and introduce the metric Cab = diag(l, —1) = c"''' to lower and raise a, b 
indices. 

Let us now proceed to the evaluation of the effective action FcTpi^a] up to the one loop 
order, which corresponds to the first order expansion of W^[Ja] in powers of h. As usual 



we solve ( p.l|) by the steepest descent method, let (p^aK^) be the solutions of the classical 
field equations, and expand the exponent in (|2.1|) about these background fields up to the 
second derivative of iS[0(j]. The integration in ( p.l| ) is now Gaussian and we can write to 
this one loop order, 

^iWiM ^ (det A,,(a;,y))-^/' (2.5) 
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where W^'^^ [Ja] = '5[</)(°)] + / dV-xJ°-(j)^^\ and Aab{x, y) is a 2 x 2 matrix defined by A+_(a;, y) 
y4_4.(a;, y) = 0, and 



^++ix,y) = 



6(t)+{x)6(f)+{y) 



AO) 



A__{x,y) 



S(j)-{x)6(f)-{y) 



AO) 



(2.6) 



In terms of tfie propagator G{x, y) = A ^{x, y) wliicli is a functional of the background fields 
ix) we can write (p. 51) as 



iy[JJ~iyW[JJ-^Tr(lnG). 



(2.7) 



The effective action which is a functional of 0^, can now be explicitly found to the same 

and the fact that 0^ differs from 0^*^) by a term of order h one can 
S[<t>a] + / dr'xJ°'(l)a, and we have finally 



order. Using ( ^.2| ), 
show that 



rcTp[0a] ^5[0J 



-Tr(lnG). 



(2.8) 



This formalism can be extended to curved spacetimes without difficulties if the spacetime 
is globally hyperbolic . The hypersufaces of constant time are now Cauchy hypersurfaces, 
the in and out states are defined in the Cauchy hypersurfaces corresponding to the far past 
and far future respectively, and the spacetime integral must be taken now with the correct 
volume element. 

Let us now compute explicity the CTP effective action for a conformally coupled massless 
real scalar field in a nearly conformally fiat spacetime. Since the detailed calculations were 
explained in Paper I, here we give only a summary of the main results. Our spacetime 
background is a spatially fiat FLRW universe with small perturbations 



(2.9) 



where rj^y = diag(— 1, +1, ■ ■ ■ , +1), exp{2Lj{ri)) is the conformal factor, r] is the confor- 
mal time which is related to the cosmological time t hj dt = exp[uj{ri)]dri, h^y{x) is 
a symmetric tensor which represents arbitrary small metric perturbations and we have 
also introduced the nearly fiat metric g^^p which is conformally related to cjj^ip. The clas- 
sical action for a free massless conformally coupled real scalar field $(x) is given by 



and 



~g^^"^^^^^^ + ^(n)i?$2 , where ^(n) = (n - 2)/[4(n - 1)], 

R is the Ricci scalar for the metric g^^. Because of the conformal coupling we can define a 
new field 0(a;) = exp[(?T,/2 — l)u;(?7)]$(a;) and the action Sm after one integration by parts 
takes the form 



Sm[g, 
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(2.10) 



which is the action for a free massless conformally coupled real scalar field 0(a;) in a spacetime 
with metric g^j,^, i.e. a nearly fiat spacetime in our case. Although the physical field is $(x) 
the fact that it is related to the field 0(x) by a power of the conformal factor implies that 
a positive frequency mode of the field 0(x) in flat spacetime will correspond to a positive 
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frequency mode in the conformally related space. Thus non trivial quantum effects will be 
due to the breaking of conformal flatness which in this case is produced by the perturbations 
h^y{x). Let us expand the above action in terms of these perturbations 

oo 

5^[V,0]=E5^^[V,</'] (2-11) 

i=0 

where the first term {i = 0) is simply the action for the field in fiat spacetime. 

To the classical action for the matter fields Sm we have to add the action of the gravi- 
tational field Sg[gfj_iy], i.e. the Einstein-Hilbert action. In order to renormalize the effective 
action later on we need to add appropriate terms quadratic in the Riemann tensor (see Eq. 
(3.6) in Paper I). We can also expand this action as Sg[gfj,u] = Y^iZo S^^^[uj, h^^], where we 
note that Sjj^^ [to] depends on to only. 

We now introduce two scalar fields (f)+{x) and 0_(x) which coincide at some future time 
T, 04.(T, x) = 0_(T, x), and which evolve in two different geometries given by /i^^ and h~j, 
such that /i^j,(T, x) = /i^^(T, x). The CTP effective action for the gravitational and matter 
action can be written as 

TcTP = Sg[uj, /i+ ] - Sg[uj, h-,] + T^[h%, 0±] (2.12) 

where contains the quantum effects of the scalar field which one can compute using 
Note that since our action is quadratic in the fields the one loop order result 
is exact in this case. Now the propagator Gah cannot be found exactly, thus we ex- 
pand Aab in the metric perturbations as Aab = + V^l^ + V^l'' + where corre- 
sponds to the Minkowski background and V^l^ are the inhomogeneous corrections, their 
explicit values up to second order are given in Sec. Ill of Paper I. Then we can write 
Gab = Gl [6,b - V}pGl - + v}^^Glv}pG% + ■■■], where the unperturbed propa- 

gator, defined by = 6ab, is a matrix such that G° ^ = Ap, G^__ = -Ad, = -A+ 

and G^_^ = A~ , where Ap and A^ are the Feynman and Dyson propagators, respectively, 
and A='= are the Wightman functions. Substituting this into ( p.8| ) and expanding its loga- 
rithmic term we get up to second order in the metric perturbations, 

rm[h%,^±] ^ E ,0+] - 5»[/.;„0_]) - '-TrilnG'J 



--Tr 
2 



2 



, (2.13) 



where we have defined v}^^ = V^!]. and V"i*^ = — V"i!l. In Paper I we did not write the terms 
which do not depend on the 0+ field because such terms were not needed to derive the field 
equation for the mean field 0[O], see (|2.4]). 

The explicit computation of the different terms was given in Paper I, the new term 
can be easily evaluated following closely that reference. After di- 



mensional regularization [|T^] of the divergent terms, renormalizing with the action of the 
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gravitational field and substituting the field equations for 0± we get the renormalized effec- 
tive action for the gravitational field up to second order as 



+ ^ I d'xd'y \m^Jx)R-'^^-^{y) - R^{x)R-^{y)\ K+(x - y- fc) 



+ - 



4 
a 

'4 

a 



d^xd'^y [3i?;,„^(x)i?-'^^"^(i/) - R-{x)R'{y) K-(x - y; fi) 
d'^xd'^y 



?+ i 



3i?+,„,(x)i?-^'^"^(y) - R-^{x)R-{y) K{x - y), 



(2.14) 



where a = (28807r^) ^, /i is a renormalization parameter, 



K^{x - y;n) 



K(x - y) 



(p2 ^ ie) 



ip-(^-y)(2vr?)0(V)^(V), 



(2.15) 



and Sj^lg^^j^] are local terms coming from the gravitational field only (see (3.29) in Paper I): 



R{x) 
IQnG 



d^x{-g{x)Y/^ 

+2a J d^x{-g{x)) 
+a f d'x{-g{x)Y/' [i?^,,^(x)/?^'^"'3(x) - R^,{x)R>^'' (x)] u;{x). 



1/2 



M\2 



(2.16) 



Here terms such as R^uaiB refer to the metric g^u and are linear in hf^u, and terms with tilde, 
RiJLva(3, refer to the metric g^j^y. 

It is remarkable that the Tctp is explicitly gauge independent since it is given in terms 
of the curvature tensor, i.e. it is invariant under the infinitesimal coordinate change x^ — >■ 

+ C^{x), for arbitrary fields C^{x), which produce the change of h^^ by /i^^, + 2^(^^j,) in 
the metric perturbations. No gauge choice was made in the above calculation. 



III. INFLUENCE ACTION FOR WEAK INHOMOGENEITIES 

As explained in the introduction we now consider the interaction of the scalar field 
(/)(x) and the gravitational field from the point of view of a quantum open system. The 
gravitational field which is the field of interest to us will be the "system" and the quantum 
scalar field will be the "environment" . Since we wish to know the effect of the environment 
on the system we will trace out the degrees of freedom of the scalar field. 

We will now summarize Calzetta and Hu |^ arguments leading to the relation between the 
CTP effective action and the influence action in the semiclassical case. We start by writting 
the generating functional J-] for a system in which we have the gravitational field, 

here symbolically called g{x), and a scalar field (t){x). We can write the classical action for 
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these two fields as: Sg[g] + Sm[4>] + Sint[g, 4>]- Note that these terms can be easily identified 
in our case from the expressions of the matter and gravitational actions. 

If we consider the field substituted by the two fields g{x) and 0(x) in (|2.1|), we can 
write 

giiy[J+,J_] ^ J 'p\^g+^'p^g-y{Ss[g+]-S,[g~]+.J+g+-J-g-+Sip[g+,g-,T]) ^ f^^^^-^ 

where we have defined 

^iSiF[g+,g^,T] ^ f 'j^^^^-j'j-)^^_-j^i{Srn[4>+]-S^[4>-]+Si„t[g+,<l>+]-Si„t[g~,<P-]) ^32) 



In ( p.2|) it is understood that we sum over all fields 0+ and 0_ with negative and positive 
frequency modes, respectively, in the remote past which coincide at time T (remote future), 
and a similar interpretation is assumed in the path integration of ( |3.1|) . Since we are only 
interested in the gravitational field we couple external sources J+{x), J-{x) to the fields 
g{x) only and not to the scalar fields Note that we are using a symbolic notation since 

the gravitational field is tensorial and so are the external sources J{x). 

The interesting point here is that ( ^.2| ) under the above interpretation is exactly the 
infiuence functional at time T, T[g^ , g~ as defined in and thus SiF[g^ , g~ ,T] is the 
infiuence action. It is the action one has to add to the classical actions Sg[g^] — Sg[g^] to 
compute any quantum probability for the transition from an initial state at t ^ — oo to a 
final state at the future time T. It is also the essential ingredient for the evolution operator, 
from initial to final time, for the reduced density matrix in this case 0. In this context the 
fields + and — and the sign difference appear as a result of the double integration needed 
to go from a transition amplitude to a probability. 

We do not attempt to quantize the gravitational field even though the formal path 
integration in seems to indicate so. Such quantization is, of course, a highly non 

trivial problem, we do not know the measure in ( p.lD and even in the linear case we have 
to deal with the gauge freedom. But if we restrict ourselves to the classical approximation 
it is easy to find the CTP effective action TcTp[g^ ■, g~] from the generating functional 
iy[J+, J_] in ( |3.1| ) using the Legendre transformation ( |2.3| ). We just need to follow the steps 
which lead from to ( p.8|) but to zero order in ^, thus we have simply VcTp[g^ , g~] = 
^gld^] ^ ^gld^] + SiF[g~^, g~, T]. On the other hand the TcTp[g^, g~] in this equation is just 
the renormalized CTP action r^j,p[^^] found in ( p.l4| ) since that action follows from the 



path integration ( |3.2| ) to one loop order for the matter fields (including scalar sources) and 
the explicit substitution of the field equations for these quantum fields. The divergencies 
are removed by appropriate terms in the classical gravitational action 5*^ as has been shown 
in Sec. |l|; note that these counterterms are implicitly assumed in (|3.1|) to remove the 
divergencies of (|37^) . Thus we have 0, 

^CTplt] = Sf[~g^]-Sf[r] + SfAh\h-,T], (3.3) 

where we have written /i^ in S^p to emphasize that it depends on the metric perturbations 
only. Comparing this expression with (|2.14|) we can write the infiuence action as: 
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J dVyC;^^p{x)C-^''-^{y)K-{x-y-fi) 
+3a J d'xd%C;^^p{x)C-^''^^{y)K{x - y). (3.4) 

Here we have introduced the Weyl tensor C^uap^x), after using the following expressions, 
which are easily shown to be satisfied for an arbitrary function f{x — y), 

J d^xd^yfix - y)[ i?^.„/3(x)i?'^^°^(y) - 4i?^,(x)i?'^^(y) + R{x)R{y)] = 0{h 



d xd yf{x - y) 



C^,^p{x)C>''"'''{y) 
-R^.ap{x)R^'''"'{y) + 2i?^,(x)i?^^(y) - ^R{x)Riy) 



OiK.), (3.5) 



where the metrics may be different at the points x and y, i.e. one may have h^y{x) and 
K^iy)^ respectively. 

Let us now write the influence action Sfp[h^J\ in a form which is more appropriate for 
the analysis of fluctuations [Q. We first decompose Sfp into its real and imaginary parts, 
because to the quadratic order in the metric perturbations the real part gives information 
on the dissipation whereas the imaginary part is related to noise, 

SfF[h%] = TiF[h%] + ifiF[h%]. (3.6) 

For this we need to decompose the kernels K^(x — y; Ji) and K(x — y) of (|2.15|) into their 
real and imaginary parts also 

K^{x -y;Ji) = 7g(x -y;Ji)± i%{x - y), 

K(a; - y) = 7o(x - y) - Z7e(a; - y), (3.7) 

where 

le{x -y;fi) = I cosp ■ {x -y)\n—, 

/i 

cosp • (x — y)0{—p^), 

sinp ■ {x — y)9{—p'^)sgn{—p^). (3.8) 

From the above definitions it is easy to show that these kernels verify the following symmetry 
relations: 

'leix - y;Ji) = %{y - x;fi), %{x - y) = %{y - x), %{x - y) = -^oiv - x). (3.9) 
It is also convenient to introduce two new kernels ^{x — y; Ji) and H(x — y; Ji) as 





r d'^p 


2 


J (27r) 




d'^p 









d'^p 


u 
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7(a; -y]fi)= %{x - y; fi) - jo{x - y)sgn{x^ - y^) = ^{y -x]fi), 

H(x -y;fi) = 7e(a; -y;Ji)+ 7o(x - y). (3.10) 

With these definitions the real and imaginary parts of ( |3.6| ) may be written, respectively, as 
3a 
Y 



d^xd^y 



nix -y- ft) + 



(3.11) 



3a 



riF[h%] = — d'xd'y C+,^(x)-C;,,^(x) U^-y) C^^'^^^iy) - C-^-^^^iy) , (3.12) 



Note that the integrand of the imaginary part depends (quadratically) only on the differences 
of the field at the points x and y as one expects of a term which is going to be the origin 
of noise. In order to write the infiuence functional in the standard Feynman- Vernon form 
0] it is convenient to change the limits in the time integrations of the above expressions so 
that x° > y^ always in the dissipation and fiuctuation terms. For this we use the following 
further symmetries of the above kernels 



y; At) = 7e(z/° - a;°,x - y;/i), 7e(x° 



7(x° - ?/°, X - y; /x) = 7(1/° - x - y; fi), 



7o(x°-l/°,x 



y) = 7e(y°-a;°,x-y), 
y) = -7o(y°-x°,x-y), 

(3.13) 

and that, for an arbitrary function f{x'^,y^), the integral J!^oo I^oo f i^^ ^ V^) 
2 I-oo dy^fix^, y^) if /(x°, y^) = f{y^, x^) and is zero if f {x^, y°) = -f{y°, x°). Since 

the kernels 7o(x — y) and 7e(a; — y) will be related respectively to dissipation and noise we 
will introduce the dissipation and noise kernels as follows 

D(a; -y) = -3a7o(x - y), 

N{x-y) = 3a%{x-y). (3.14) 
Finally we can write the infiuence functional in the standard form as 

S^fK] = ^ / d'xd'y^ix - y; ^) fc+ (x)C+'^^"^(i/) - C-Jx)C-^-^{y) 



dx" / dy" / dVyAC,,^p{x)D{x - ynC^^'^^iy)} 

/OO PX*^ p 

dx'> / dy° / dVyAC^,o.p{xMx - i/)AC"^^"^(i/). (3.15) 
-00 J—oo J 

where we have used the notation 

AC{x) =C+{x) -C-{x), {C{x)} = C+{x) + C"{x), (3.16) 

for the difference and sum of the + and — fields at the same spacetime point. 

The first term in the infiuence functional is non-local and will contribute to give a non 
local stress tensor in the semiclassical equations but does not mix the + and — fields. On 
the other hand, in the second and third terms, which are also non local, there is mixing of 
the + and — fields. The second will contribute to the dissipation and the third to noise as 
we will see in the next section. We should remark also the explicit gauge invariant form of 
the infiuence functional above. 
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A. The fluctuation-dissipation relation 



It is easy to derive a relation between the dissipation and noise kernels in our case in 
analogy with the quantum Brownian model of ref. [0 or the case of a quantum scalar field 
in an anisotropic Bianchi I cosmology [Q. 

We first note that the dissipation kernel D(x) of ( p. 14 ) can be written as a time derivative: 



D(x) = ^7(^), (3.17) 

where 

where l/|p°| in this integral must be understood as its Hadamard's finite part distribution. 
Then the fluctuation-dissipation relation takes the form 

N(a;) = J d^x'KFDix-x')-f{x'), (3.19) 

where the fluctuation-dissipation kernel Kfd{x — x') is given by the distribution 

KFDix-x') = 5%x-x') -^q" COS q^iri-r]'). (3.20) 

JO vr 

To check this fluctuation-dissipation relation one may simply substitute ( p.2(]| ) into the right 
hand side of ( |3.19|) , perform the integrations and the result is the noise kernel deflned by 



]8[) and ( |3.14| ). This relation is important since it gives a direct connection between the 
effect of quantum fluctuations of the environment on the dissipation of the gravitational 
fleld inhomogeneities. 



IV. EFFECTIVE INFLUENCE ACTION AND STOCHASTIC SEMICLASSICAL 

EQUATIONS 



In this section we obtain from the influence action of Sec. |T| the effective influence 
action. From that effective action the explicit semiclassical Einstein-Langevin equations will 
follow. 

We start by recalling some well known relations of frequent use in statistical 
physics. The path integral Gaussian identity [Q, / X'^exp[— 1(^, L^) + (&, + c] = 
(detL)-i/2gxp[i(6,L-i6) - c], where L is a linear operator acting on the fleld ^(x) and 
the brackets in the exponents stand for (^, C) = / ^{x)({x)d"'x, can be used to show that 
(change L by A~^, b by ik and take c = 0) 

^k] = e-^C^.^'^) = f V^V[^]e'^^'^\ (4.1) 
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where V[$] is given by 



If we interpret V[^] as the probabihty distribution functional of the field ^{x) then $[/c], 
which is its functional Fourier transform, is the characteristic functional. The mean value 
of a given functional A[^] is defined by 

(A)^ = I V^rmi^] (4.3) 

and the n-correlation functions of the field ^(x) may be derived from the characteristic 
functional as 



(4.4) 



k=0 



Since in this case the characteristic functional is Gaussian we have that {C,{x))^ = and that 
the two point correlation function is {^{x)C,{x'))^ = A{x,x'). 



Let us go back now to the influence functional derived in Sec. Ill, 



TiF[h%] = e^^^^t'^'"'] = e^^ipif^'^-^^e'^'^^^^^^ (4.5) 

where the real and imaginary part of the influence action Tjp and Tjp, respectively, are 
given in ( |3.15| ). We recall now the observation that the imaginary part of the influence 
action depends on the difference between the + and — fields only, more precisely, it can 
be written as f7F[/i^^(x)] = tiF[AC'"'''^{x)]. This is the signal g] that the effect of the 
environment on the system, given by this part of the influence functional, is equivalent to a 
classical stochastic source ^^ua/si^) (a tensor field in this case) whose probability distribution 
has such part of the influence functional as its characteristic functional. In fact, using ( |4.1|) 
we see that (14.51) can be written as 



^iF[h%] = J V^V[^]expt{tjp[h^,] + J d'x^,,^^{x)AC''^-^{x)Y (4-6) 
where V[C] is the Gaussian probability distribution given by (see (^]2|) and ( p.l5| )) 

-i /dWs,C(x)(N(x-j/))-iC(j/) 

Therefore using ( [4.3| ) we can interprete the influence functional as the following mean value 

^iF[hU = (e^^- (4.8) 
where the effective influence action is defined by 

Sf/[h%, e] = TifK] + J d'x^,^^p{x)AC^''-P{x). (4.9) 
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The effect of the environment (quantum fields) on the system (the gravitational field) is 
completely characterized by this effective action, the tensor i^j.yap{x) plays the role of a 
stochastic source with the Gaussian probability distribution given by (|4.7|). This tensor 
has the symmetries of the Weyl tensor, i.e. it has the symmetries of the Riemann tensor 
and vanishing trace in all its indices. The kernel N, which appears in ([4.7| ), can thus be 
interpreted as the noise kernel in our problem. Since the probability distribution is Gaussian 
the noise kernel is the two point correlation function of the stochastic source. This source, 
in fact, is completely characterized by the relations 

{i^,ucp{x))^ = 0, 

{^t^uaf3ix)^p^xeiy))^ = T^i.a/3paAeN(x - y), (4.10) 

where the explicit form of the tensor T^ua/Spaxe is given in the Appendix, it is the product of 
four metric tensors, in such a combination that the right-hand side of the equation satisfy 
the Weyl symmetries of the two stochastic fields on the left-hand side. It is easy to obtain 



these relations using ( |4.4| ), note that the characteristic functional has the form $[A;^!/o/3(a;)], 
where kpuap{x) has the symmetries of the Weyl tensor. It should be now clear also that D 
is the dissipation kernel for this problem since it is related to the noise by the fluctuation- 
dissipation relation (|3.19|) . It is worth to mention that the association of the imaginary terms 
of the effective action (or the influence action) as terms coming from the interaction of the 
field with stochastic sources has been used also in the context of scalar fields in interaction 
with other fields, or selfinteracting, in order to study the nonequilibrium dynamics of these 



quantum fields [IS . 



A. Stochastic semiclassical equations 



We are now in the position to derive the semiclassical equations for the gravitational 
field due to a quantum scalar field. We recall that the effect of the quantum field on the 
gravitational field is given by the effective influence action ( [4.9|) . Therefore the total effective 
action, which includes the action of the gravitational field plus the previous effective influence 
action is given by 



(4.11) 



where S^[gfj_y\ is the renormalized action of the gravitational field (|2.16|) . Since this is 
an effective action for the metric perturbations the field equations for the metric can be 
derived in a similar way as (|2.4|) . It is usually convenient to introduce new variables h^^j = 



{h^u + h/iiy) /2, A/i^i^ = /i+j, — /i^j,, the average field and the difference field respectively. 



eff 



Then the field equations are 
form, which is of more practical use to us 



0; or in an equivalent 



6h 



flU 



0. 



(4.12) 
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The functional derivations needed in the first two terms can all be found in Appendix E of 
Paper I. For the new term, we can easily prove that 



/ 



OLUTp 



X 



-2 1 d^xd^d,e''^p{x)hM 



(4.13) 



where in this expression we have assumed that the tensor field C,ij.uai3 has the symmetries of 
the Weyl tensor. 

Finally the semiclassical equation ( [4.12 ) can be written as 



1 



3a 



12 



^(0) 



^(1) 



a 
2 



+- 



'(1) 



(4.14) 



where the (0) and (1) subindices mean the zero and one orders, respectively, in terms of the 
perturbation /i^j,. Terms with and without tilde refer to tensors obtained with metrics g^y 
and g^vi respectively. G^'^{x) is the Einstein tensor, B'^'^{x), A'^'^{x) and H^^{x) are given 
in Paper I as: 

5^'^(x) = ^^''^i?^ - 2RR>"' + 2i?''^^ - 2g>"'UgR^ 



-RR 



2UgR^^ + '^-R^^^ + ^-g'^'UgR, 



H^^ix) = -R^'^Ra," + '^RR^" + ^g^^'RapR"^ - ^g^^'R^- 



(4.15) 



The tensor F^^{x) results from the functional variation of the stochastic term, it can be 
read off directly from ( |4.13| ), 



F>"'{x) 



XI 



(4.16) 



and it is symmetric and traceless, i.e. F^'^(x) = F'^^{x) and F^^{x) = 0, which means that 
there is no stochastic correction to the trace anomaly. 

From equations (|4.14| ) we may now define the effective stress tensor T^ff, which is the 
term to write on the right hand side of Einstein equations, i.e. let us write (|4.14|) as 



(4.17) 



where {T^'')q is the (quantum) vacuum expectation value of the stress tensor of the quantum 
field up to first order in h^y and we have added a classical stress tensor T^^ to include the 
case in which there is also a classical source (this was not considered from the beginning for 
simplicity). The quantum stress tensor is given by 
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(7^(1)). 



a 



a 



^(0) 



6^(0) 



{Hi;'' 



(1) 



(4.18) 



This tensor was already given in Paper I and was first computed by other means in Refs. 

. Now we have derived a stochastic correction to this tensor which accounts for the noise 
due to the fluctuations of the quantum field. 

If we now take the mean value of equation ( [4. 171) with respect to the stochastic source ^ 



we find that, as a consequence of (|4.1CI| ) 



(4.19) 



and we recover the semiclassical Einstein equations of Paper I. 

The stochastic correction to the stress tensor has vanishing divergence to the first order 
in the metric perturbations. In fact, using that g^i, = e^^g^v and that F^'^ is symmetric 



and traceless, it is easy to see that (e ^'^F'^^) = e ^"^VyF^" , then from ( [4.16|) and the 



symmetries of ^ we obtain that V yF' 



0{h 



It is thus consistent to write this term 
on the right-hand side of Einstein equations and consider it as a correction of order higher 
than (Tjp, (note that V.(T5P, = 0{hl,)). 

To summarize, the stochastic semiclassical equations ([4.17|) can be called the semiclassi- 
cal Einstein-Langevin equations for weakly inhomogeneous spatially flat cosmologies in the 
presence of a conformally coupled massless scalar field. The usual semiclassical equations 
can be seen as the mean value of these equations with respect to a tensor field source ^ 
with Gaussian probability distributions ([4.7|). This stochastic source couples to the Weyl 
(conformal) tensor of the spacetime metric in the form given by ( |4.9| ), this means that it 
has the symmetries of that tensor and thus that it has only 10 independent components at 
each spacetime point. 

The fact that the stochastic source couples to the conformal tensor should not come as 
a surprise since we expect that, for a conformal quantum field, non trivial quantum effects 
should be a consequence of breaking the conformal symmetry of the spacetime, which is 
characterized by the conformal tensor. For instance, it is known that the probability density 
of pair creation in this case, or in the presence of only small anisotropics, is determined by 
the square of the Weyl tensor Thus as it has been emphasized in there is a direct 
relation between particle creation and noise. 



V. METRIC FLUCTUATIONS IN FLAT SPACETIME 

In this section we make a simple application of the semiclassical Einstein-Langevin equa- 
tions obtained in the previous section to the case in which the background spacetime is not 
cosmological, i.e. when u = in ( |2.9| ). This restriction simplifies considerably the semi- 
classical equations. We will take here a perturbative approach in which the semiclassical 
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corrections to Einstein's equations are seen as analytic perturbations (in h) to the classical 



Einstein's equations; see Simon |^ for a justification of this point of view. 

Let us assume that we have some weak stress tensor source T^^ in fiat spacetime. For 
instance, we could have a cosmic string or a bulk of Newtonian matter. Such a source will 
produce some classical linear inhomogeneities /i^^, and the spacetime metric will be 

gt^u = Vi^u + hl^u- (5-1) 

If we have a massless conformally coupled quantum scalar field in this background, the stress 
tensor acting on the spacetime will now be T^^, + T^^-'^. The semiclassical equations up to 
first order in the metric perturbations can be written as 



rM+T^"{x)], (5.2) 



where T^;^-^ follows from ([4.17|) and ( [4.18|) with u = 0; note that in this case (T^op? ~ ^■ 
The stochastic term 2^^*^ will produce a stochastic contribution /i^*^, to the spacetime 

inhomogeneity, i.e. h^y = h^^^ + h'^^^. Let us now consider such contribution. 

Substituting into ( |5.2| ) and taking into account that /i^^ already satisfies the classical 

equation, we obtain a linear equation for the stochastic term h^^^, which we may write as 

S^t = 2F^"' = -Ad^dpe""^, (5.3) 

together with the harmonic gauge condition, (/i**^ — ^rj^^h'^^y^ = 0, which was used to write 
the previous equation. The solution of these equations as a Cauchy problem with boundary 
conditions /;.^*^(— oo, x) = 9j/;,^*^(— oo, x) = 0, is given by 



/ij,(t, x) = WnG^ f dt' [ d'^Dnit - t', x - x')^^t(^', x'), (5.4) 



dt' , 

where Dji is the retarded Green's function 



Dnix-x') = -—^ -(5(t-t'- |x-x'|). (5.5) 

47r|x — X I 

Note that after imposing the above boundary conditions for /i^*^, the gauge has been com- 
pletely fixed because two metric components within the harmonic gauge can differ by 2(^^^i^ 
where is an harmonic vector field but this vector field is zero when such boundary condi- 
tions are imposed. Since S"^^ is linear in the stochastic source, it is obvious from (|5.4| ) that 

Let us now compute the two point correlation function of the stochastic metric fiuctua- 

,(g-t(t-|x-x-|,xOgf,(s-|y-/|,yO)^ 



tions h^^ 



Jr3 |x — x'lly — y' 



(5.6) 
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here we made use of ( ^.41) and performed two time integrations. Now since 5**^ has the form 
of a hnear operator acting on the tensorial stochastic source, ^, the correlation function of 
S"*^ can be written in terms of the correlation function of that source, i.e. in terms of the 
noise kernel N. After some simple manipulations we obtain 

{S7{^)s^^y)h = ^Q""'' I (^^"■^^( V), (5.7) 

where we have defined = — y'^, and Q^^^^ is the operator 

with P'^'^ = rj^^n — d^d'^, where all derivatives are with respect to a^. Substituting this into 
(^ leads to 

{h^;^it,^)hl',{s,y))^ = 67r'a{32G^fQ,^,e J ^^( V)e^^'^A+(p)A-(p) (5.9) 

where are the Wightman functions. Finally, after some further manipulations the cor- 
relation function of /i^t(^) can be written as: 

(^^t(t,x)/ife(s,y))g = {lenG^fe^^xe J f/'^K,,,(a - z)N(z), (5.10) 

where Kreg{x) is the appropriately regularized kernel defined by the Fourier transform in- 
tegral / d'^pexp^ip ■ a;)A+(p)A~(p) = K^egix). Since the noise kernel is non-local, see 
equations ( p.8| ) and ( |3.14| ), the noise is colored. A simple order of magnitude estimate 
of the above result gives that it is of the order of (Planck time/time interval)^, as one 
should expect of quantum fluctuations. Thus if we measure the gravitational potential at 
a given point at different time intervals we should find variations in this field of the order 
v(^oo^oo) ~ (Planck time/time interval)^. 



APPENDIX A: 



The T^uafBpaxe tensor is given by 



+4 

-3 



'np[^iVu]aVX[aVl3]e + Vp[aVl3]<7VX[pVu]e + Tla[piTlv]l3VX[pVa]i 

Vp[pVl3]ailx[ar]u]e + rjp[pVa]arix[urili]e + Vp[vV<J\arix[/3'np]e + VpWVulaVXiaVple 

Vlia {TlpXTlcT{vTlf3)e + 'n^eVp{iyVl3)X - 'naXVp{iyVl3)e - VpeVa(iy1lf3)x) 
{jlpXVcT{pTla)e + V<Ter]p{f,r]a)X - TlcTXTlp{t^Va)e - TlpeV<T{iiTla)x) 

~Vua {ripxria{iirip)e + r]„eVp{pVi3)x - V<TxVp{pVi3)e - r]peV'7{iJ.Vi3)x) 

)]}' (Al) 

it comes from the functional derivative of the Weyl tensor with respect to itself taking into 
account all its symmetries. 
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